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Abstract. We construct a seven-component model of the in-host dy¬ 
namics of the Human Immunodeficiency Virus Type-1 (i.e, HIV) that 
accounts for latent infection and the propensity of viral mutation. A 
dynamical analysis is conducted and a theorem is presented which char¬ 
acterizes the long time behavior of the model. Finally, we study the 
effects of an antiretroviral drug and treatment implications. 

1. Introduction and Model. The human immune system is a complex 
network of interacting cells, cell products, and cell-forming tissues that pro¬ 
tects the body from pathogens, destroys infected and malignant cells, and 
removes cellular debris [11]. A key component of this system is the white 
blood cell known as a T cell. In particular, the CDd"*" T cell moves through¬ 
out the body, identifying bacteria and viruses and directing the immune sys¬ 
tem’s attack. Within healthy individuals the blood concentration of these 
cells is relatively constant, around 10® cells/ml [5]. HIV molecularly recog¬ 
nizes CD4'*'T cells as compatible for viral replication. As HIV is a retrovirus, 
it replicates within a host via reverse transcription using its RNA and an 
enzyme (reverse transcriptase) to create a strand of HIV DNA, called a 
provirus, which carries its genetic information. Once the provirus is created 
within a CD4’*'T cell, its immunological function ceases and it is used to cre¬ 
ate new virions [1]. Interestingly, reverse transcription is quite error prone, 
and the probability of mutation is high (around 22% of infected cells carry 
proviral genomes with at least one mutation [8]). Without being activated 
by antigens, an infected CDd+T cell, containing a provirus, fails to produce 
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new virus particles. In such a case, the T cell is referred to as latently 
infected [2, 9, 10]. 

The following model describes interactions between T cells and the popu¬ 
lations of a wild-type and mutation of HIV. It incorporates latent infection 
of T cells and accounts for a single mutation within the virus. Specifically, 
the model is comprised of seven nonlinear ODEs given by 


dt 

dis 

dt 

dl/g 

dt 

dK 

dt 

dir 

dt 

dLr 

dt 

dVr 

dt 


X - drT - ksTVs - KTVr 


(1 - p)(l - fi)ksTVs + asLs - djls 


p{l - p)ksTVs - asLs - dLLs 


Ngdils — dyVs 




(1 — p)lJ,ksTVs -t (1 — p)krTVr + arLr — dflr 


PfiksTVs + pkrTVr — arLr — dLLr 


Nrdjir — dyVr- 


(7CM) 


Here T denotes the population of healthy T cells whose maturity rate is A 
and death rate is dp- The parameter kg is the rate of infection by the wild- 
type virus, whose population is denoted by Vg. Note that the T cells infected 
by the original strain, kgTVg, contribute to several other populations. Ig and 
Lg denote the populations of wild-type actively or latently infected T cells, 
respectively. The proportion of T cells which become latently infected is 
given by p. Ir and Lr are the populations of T cells that are actively and 
latently infected by the mutated strain, while the mutation rate of the wild- 
type virus is /i G [0,1]. Note that T cells which are infected by the wild-type 
strain may create a mutated provirus and thus the terms (1 — p)iikgTVs 
and ppkgTVg are included. The activation rates of latently infected cells to 
actively infected are ag and ar for the wild-type and mutated strains. The 
death rates of actively and latently infected T cells, regardless of strain, are 
dj and dp ■ Ng and Nr are the numbers of virions (wild and mutated) released 
by an infected T cell over its lifespan. Finally, dy is the viral clearance 
rate for both strains. See the review [6] for more information. Parameter 
values are in Table 1, while Fig. 1 presents a simulation of (7CM) with this 
data. Both strains persist in the simulation, but the wild-type dominates 
the mutated strain by five orders of magnitude. In subsequent sections we 
will utilize analytic and computational methods to study (7CM) and discuss 
the effects of antiretroviral drugs on these populations. 


2. Analysis and Dynamics. With the model finalized, we turn to an anal¬ 
ysis of its behavior. First, well-posedness of solutions to (7CM) can be shown 
using the Picard-Lindelof Theorem and known a priori bounds. The details 
are standard so we omit them. Instead, we focus on the asymptotic behavior 
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Table 1. Variables and Parameters 



Variables/Parameters 

(Initial) Value 

Reference 

T 

Uninfected CDd"*" T cells 

10® mpi 

[5] 

Is 

CDd"*" T cells actively infected by wild-type virions 

0 


Ls 

CD4"'' T cells latently infected by wild-type virions 

0 


Vs 

Wild-type HIV virions 

10"® mpi 

[4] 

Ir 

CD4"'' T cells actively infected by mntated virions 

0 



CDd"*" T cells latently infected by mutated virions 

0 


Vr 

Mutated/Drug-resistant HIV population size 

0 


X 

Rate of supply of immunocompetent CDd"*" T cells 

10"^ ml~^ day”^ 

[8] 

(It 

Death rate of uninfected T cells 

0.01 day-^ 

[8] 

ks 

Infection rate of T cells by wild-type HIV 

2.4 X 10“® ml day" 

"1 [8] 

kr 

Infection rate of T cells by mutated HIV 

2.0 X 10“® ml day" 

[8] 

P 

Rate of latent infection in T cell population 

0.1 

[2] 

di 

Death rate of actively infected T cells 

1 day“^ 

[8] 

dL 

Death rate of latently infected T cells 

4 X 10"® day"^ 

[2] 

as 

Activation rate for latently infected wild-type T cells 

0.01 

[2] 

QJt’ 

Activation rate for latently infected mutated T cells 

0.01 

[2] 

Ns 

Burst size of wild-type virus 

3000 

[8] 

Nr 

Burst size of mutated/drug-resistant virus 

2000 

[8] 

P 

Mutation rate from wild-type to mutated 

3 X 10"® 

[8] 

dv 

Clearance rate of free virus 

23 day“^ 

[8] 


of (7CM) as t —7- oo. A number of long calculations show that exactly three 
steady states exist - Eq, Er and E^, given by 


En — 


1 — \ 


d'j’ 


0 


0 


0 

Er = 

0 


0 


i 0 ; 



d'j'Rj- 

0 

0 

0 


where the newly-introduced parameters Rr, Rs, Rc, c and cic are defined as 


V 


E, = 


dx Rc 

d^dyia^-l) _,N 

d I kg Ng (1—/^)(CT— 1)^ ^ ' 

dy (g-e-l) 


ax’l^O'c —Ij / r> _ 1 \ 

fcs(l-Ai)(o--l) 

_ ^td^dy _ ( ry _ 1 ^ 

dikrNr{l-lJi){cT-l) 

Rr[dL+ar)0--tJ.)(<y-l) ~ 1) 


Rr = 


\krNr{ar + (1 — P)dL) 

dTdv{dL E CXr) 


Rx = 


R. 


XksNsjas + (1 -p)dL) 

dxdvidi + tts) 


Rc = {l-p)Rs, cr=—, CTc = (l-^)iT. 

■t\'p 


Eq is the infection-free steady state, in which all populations but T are zero. 
Within Er, the mutation dominant steady state, only the mutated strain re¬ 
mains present as R, Lg and 14 are zero. Ec is the state in which both strains 
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(a) Healthy T cells (b) Wild-type HIV Strain 



(c) Mutated HIV Strain (d) Total Virus Load 


Figure 1. A log-Iog simulation of (7CM) with parameters given in 
Table 1. Notice that the wild-type viral strain is dominant. 

persist, and hence is referred to as the coexistence steady state. To repre¬ 
sent biological quantities, all steady state components must be nonnegative. 
Clearly, this is the case for Eq. Since all parameters are positive, we see 
that Rr > 1 implies that Er satisfies this requirement. Similarly, Ec requires 
Rc > I and cic > 1, which implies Rg > Rg > Note that 

(Tc > 1 implies a > 1. Fig. 2 displays the regions in the -Rs)-plane for 
which steady states are nonnegative. 

Finally, we wish to determine whether such states correspond to locally 
asymptotically stable (l.a.s) or unstable equilibria. In this vein, we prove 

Theorem 1. The stability properties of the equilibria depend only upon the 
two parameters Rr and Rg. Moreover, 

1. The infection-free steady state Eq is l.a.s. if and only if 

Rg<l/{l — p.) and Rr < 1, 
and it is unstable otherwise. 

2. The steady state with only the mutated virus, Er is l.a.s. if and only if 

Rg < Rr /(1 — p) and Rr > 1, 
and it is unstable otherwise. 

3. The coexistence steady state Ec is l.a.s. if and only if 

Rg>l/{l — p) and Rs>Rr/{l — p), 
and it is unstable otherwise. 
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Rs 



Figure 2. Regions of biologically-relevant existence for steady states 
in the (RrjRs)-plane. 

Proof. As in [3, 7] we consider the linearization of the system to justify 
conclusions about the its local asymptotic behavior. To begin, we consider 
the state Eq, and evaluate the Jacobian of the system at this point. In 


particular, we 

are 

interested in 

the eigenvalues 

of the 

matrix 




/ 

-dj^ 

0 

0 

d'jp 

0 

0 

Xkr 

d'j’ 

\ 



0 

-di 

as 

d'j' 

0 

0 

0 




0 

0 

—dL — 


0 

0 

0 


Vf(Eo) = 


0 

diNs 

0 

— dy 

0 

0 

0 




0 

0 

0 

{p—l)Xfj,ks 

d'j’ 

-di 

ar 

{p—\)Xkr 

d'j’ 




0 

0 

0 

pXpks 

djp 

0 

— dL — O-r 

pXkj^ 

d'j' 



1 

0 

0 

0 

0 

dj Nr 

0 

—dy. 

/ 


The seven eigenvalues are Ai = —(It', ^2,3,4 defined by the roots of the cubic 

030 ;^ -|- a2x‘^ + aix + ao 

where 03 = 1 and 02 = djdr + didT + Urdr -h dxdy are both positive, and 
oi = dT'^idy -|- dj){dL ctr) dx'^dvd] — dTdj{l — p)XkrNj., 

CLq = dT^d\/d/(dT -h (Xr) — djdT'^XkrNr{(^r T (1 ~ P')dL) 
and A5^6,7 defined as the roots of the cubic polynomial 

h'iX^ + b2x‘^ + bix + bo 

where 63 = 1 and 62 = djdr + didx + ocgdr + drdy are both positive, and 
61 = dT‘^{dy + dj){dT -\- cxs) + dT^dydi — dTdj{l — /i)(l — p)XksNs), 
bo = dr^dydiidi + ««) - - iJ.)XksNs{as + {1 - p)dL)- 
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(a) Behavior in Region 1 





HeatfilyT-cells 



Figure 3. stability of steady states in various regions of the (RrjRs)- 
plane. In each figure, the T component of Eq, Er, and E^ is represented 
by dotted, dashed, and dot-dashed lines, respectively, while the solid 
curve is the solution T{t) with T(0) = 10®. 

It follows from the Hartman-Grobman Theorem that Eq is l.a.s if and only 
if Re{\i) < 0 for all eigenvalues, Aj, of Vf(Fo). By the Routh-Hurwitz 
Theorem, the roots of the cubics possess negative real parts if and only if 
“0; “1,02,03 > 0 and 0201 > 0300- Clearly, 02,03,62,^3 > 0 . For the first 
cubic, oi,ao > 0 and 0201 > oq are satisfied for < 1. For the second 
cubic, the corresponding requirements are satisfied if and only if Rg < 
Hence, Eq is l.a.s for the specified conditions. A similar analysis proves the 
second and third conclusions of the theorem. For brevity, we omit the details 
but take a computational approach to demonstrate their validity. 

Fig. 2 displays the regions in which nonnegative steady states exist. In 
the area of the (R,., i?s)-plane where Rg > and R,- < 1 (Region 1), 
both Eq and Ec exist. For Rg > y^Rr and R^ > 1 (Region 2), all three 
steady states exist. For Rg < y^Rr and R^ > 1 (Region 3), both Eq and 
Er exist. Simulations are used to demonstrate the stability of coexisting 
states. We first begin with parameter values placing {Rr,Rs) in Region I. 
Note in Fig. 3(a) that T{t) (solid) tends to the corresponding value of Ec 
(dot-dashed). Next, (R^, Rg) is adjusted to lie within Region 2, in which all 
steady states exist. Fig. 3(b) shows that T{t) still tends to the corresponding 
value in Ec- Thus, simulations validate the results of the theorem, namely 
that within Regions 1 and 2, Ec is l.a.s. while Eq and Er are unstable. 
Using {Rr,Rs) within Region 3, Fig. 3(c) demonstrates that T{t) tends 
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Figure 4. Regions of stability of equilibrium points in the (-Rr,Rs)-plane 

to the corresponding value of Er (dashed) as t —)• oo. Thus in Region 3, 
Ej. is l.a.s. While only the healthy T cell population is shown, all other 
populations display the corresponding long term behavior. These findings 
and the theorem are summarized by Fig. 4. 

□ 

3. Drug Therapy. In the previous section it was shown that the long term 
behavior of (7CM) depends crucially upon parameter values, and even in 
the presence of viral mutation, the virus may be completely cleared from 
the body. Contrastingly, the mutated strain may dominate the dynam¬ 
ics, which can lead to drug resistance. Next, we consider the introduction 
of antiretroviral therapy (ART), specifically reverse transcriptase inhibitors 
(RTIs). These drugs hinder the replication process and decrease the rate of 
infection, thus altering kg and ky. Let Eg represent the efficacy rate of the 
RTI on the wild strain. We assume the RTI combats this strain, while the 
mutated strain is more resistant to the drug. Hence, the efficacy of the RTI 
on the drug resistant strain is given by Sy = aSg, where a G (0,1) represents 
the resistance level. Note that smaller values of a correspond to greater 
levels of resistance. Including ART within the previous model yields: 

= A — dj'T — (1 — Zs)ksTVs — (1 — €r)krTVr 

= (1 - p){l - p)(l - es)ksTVs + oisLs - djis 

= p(l — p){l — €s)ksTVs — OisLs — dLLs 

= Nsdjis — dyVs 

= (1 - p)fl{l — es)ksTVs + (1 - p)(l — er)krTVr + OirLr — djir 
= p/i(l — €s)ksTVs + p{l — £r)krTVr — arLr — dj^Lr 
T'^'f'djl'p — dy\^j'. 


dT 

dt 

dh 

dt 

dLs 

dt 

dVs 

dt 

dif 

dt 

dLf 

dt 

dVr 

dt 


> 


(7CM,) 
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(a) a = 0.5 


(b) a = 0.2 


Figure 5. For visual clarification the difference between R% and 7?^ 
has been enlarged by taking /r = 3 x 10“^ in Note the threshold 
differences between a = 0.5 (left) and a = 0.2 (right). 


Based on similarities between (7CM) and (7CMe), the steady states are 
merely EQ,Er, and Ec under the transformations ks (1 — es)ks and kr e-)• 
(1 — £r)kr. We denote them by Eq, Ef., and E^, respectively. The new 
parameters and (as well as R^, and a^) are defined similarly: 

Rf. = {1- £r)Rr, Rl = (1- £s)Rs- 

Hence, the stability properties can be deduced from Theorem 1 by substi¬ 
tuting Rr = R^, and Rg = Rl- 

We examine how the reproductive ratios are affected by ART for increasing 
drug efficacy rates assuming a particular resistance level, in this case a = 
0.5. In Fig. 5(a), the efficacy rate £s is plotted against reproductive ratios, 
Wg and R^, along with R^ and cr^ to view changes in long term behavior. 
The threshold values constituting these changes are denoted by £i and £ 2 - 
For values of £g < £ 1 , we observe R^ > 1 Rg > and > 1 
Rl > j^Rr- Thus by Theorem 1, we conclude that solutions tend to E^. 
Contrastingly, for ei < < £2) ^ cr^ < 1 thus solutions tend to 

E^. It is not until £g > £ 2 , when < 1 and < 1, that the efficacy is 
large enough to force viral clearance. 

To investigate the change in the long term behavior of the viral loads, we 
graph the steady state values of W and Vg separately against £s. Figs. 6(a) 
and 6(b) display the change in the drug resistant virus population. Fig. 6(a) 
shows the dramatic increase in the drug resistant population near £g = £ 1 , 
but this scale cannot accurately display values for £g < £ 1 . Thus Fig. 6(b), 
an enlargement, shows the slow increase in this population. In Fig. 6(c), 
we note the steady decline in the wild population as £<j increases and its 
extinction at £* = £ 1 . Fig. 6(d) displays the total virus population. 

The question arises whether a value of a exists which prohibits viral clear¬ 
ance altogether. A higher resistance level of the mutated strain is chosen, 
and the steady state behavior is examined. Fig. 5(b) shows the reproduction 
numbers assuming a = 0.2. The threshold £1 = 0.5 when a‘^ = 1 indicates 
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(a) Vr steady state (b) Vr steady state, £« < ei 



(c) Va steady state (d) Total virus steady state 

Figure 6. steady state values for both wild type and drug resistant 
virus population with respect to £a, assuming a = 0.5. 


the shift from the coexistent state to drug resistant dominance. There is no 
longer an 62 threshold, which implies that no drug efficacy could be large 
enough to counter the resistance level of the mutant strain. Thus, the RTI 
would clear only the wild type virus. Fig. 7(a) demonstrates that even when 
Ss = 1, the resistant steady state remains prohibitive. In Fig. 7(c), we note 
the faster decline in the wild population as increases and the earlier ex¬ 
tinction at El as compared to Fig. 6 (c). Fig. 7(d) displays the total virus 
population, demonstrating continued viral persistence as e* —1, unlike Fig. 
6 (d). Therefore, viral mutation and drug resistance can completely inhibit 
the treatment ability of RTIs. 


4. Conclusions. Upon formulating a new in-host model of HIV dynamics 
that incorporates latent infection and viral mutation, the local asymptotic 
behavior of the model was characterized and validated using computational 
means. The effects of RTIs were elucidated under differing drug resistance 
levels for the mutated strain. In particular, our analysis determined that a 
sufficiently high resistance level arising from viral mutation could render a 
persistent infection untreatable. This implies that RTIs alone are not strong 
enough to counteract the emergence of drug resistant strains. Additionally, 
as HIV possesses a strong propensity to mutate, due to common errors in 
reverse transcription, the effects of viral mutation are significant within the 
dynamics of the virus. 
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(a) Vr steady state 


(b) Vr steady state, Ss < £i 




(c) Vs steady state 


(d) Total virus steady state 


Figure 7. Dynamics of steady states for wild type and drug resistant 
virions with respect to drug efficacy, assuming a = 0.2. 
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